Suppression of weed and insect populations by living and straw mulches in sesame (Sesamum indicum L.)

In order to evaluate the effect of different weed management treatments on weeds, pest and natural enemies populations in sesame (Sesamum indicum L.), a 2-year study was conducted in East Azarbaijan, Iran in 2020–2021. The study was conducted based on randomized complete block design with four replications. The weed management treatments consisted of trifluralin use (960 g ai ha−1), wheat straw mulch (WSM), living mulches of fenugreek (Trigonella foenum-graecum L.) (FLM), bitter vetch (Vicia ervilia L.) (VLM), calendula (Calendula officinalis L.) (CLM) and one-time hand weeding (OHW). The effect of weed management treatment was significant on densities of insect pests, natural enemies and weed and also weed biomass and sesame seed yield. The lowest densities of insect pests including Myzus persicae, Brevicoryne brassicae, Helicoverpa armigera and Spodoptera exigua were observed in CLM treatment. Also, the highest densities of natural enemies Coccinella septompunctata, Coccinella undecimpunctata and Orius niger were observed in CLM treatment. The highest reductions in grass (51.0%), broadleaf (72.0%), and total (62.6%) weed biomasses and highest seed yield (1456 kg ha−1) were obtained in OHW. The seed yields in CLM and WSM treatments were not significantly different with trifluralin treatment and could be recommended in sustainable production of sesame.


Population of insect pests
Myzus persicae density was affected significantly by year (p ≤ 0.05) and weed management (p ≤ 0.05) (Table 1).The Myzus persicae density in 2020 was higher than 2021 (Table 2).The highest density of Myzus persicae was observed in weed-free treatment and decreased significantly at all treatments and the lowest value was observed in CLM treatment (Table 3).The density of Myzus persicae in FLM, VLM and OHW treatments were not significantly different.The Brevicoryne brassicae density was affected significantly by year (p ≤ 0.05) and weed management (p ≤ 0.01) (Table 1).Brevicoryne brassicae density in 2021 decreased significantly compared with 2020 (Table 2).The highest and lowest densities of Brevicoryne brassicae were observed in weed-free and CLM treatments, respectively (Table 3).The densities of Brevicoryne brassicae in trifluralin and weed-infested treatments, were not significantly different.The Brevicoryne brassicae density in WSM treatment was higher than FLM, VLM and OHW treatments.
The effect of year was not significant on densities of Helicoverpa armigera and Spodoptera exigua (Table 1).The effect of weed management treatment was significant (p ≤ 0.01) on densities of Helicoverpa armigera and Spodoptera exigua.The highest densities of these insect pests were observed in weed-free treatment (19.3 and 28.9  no.plant −1 , respectively) and the lowest values (5.4 and 6.2 no.plant −1 , respectively) in CLM treatment (Table 3).
Table 1.Analysis of variance for effect of weed management treatments on densities of insect pests and sesame seed yield.ns, * and **: non-significant and significant at p ≤ 0.05 and p ≤ 0.01, respectively.There were no significant differences in densities of Helicoverpa armigera and Spodoptera exigua in VLM, FLM and trifluralin treatments.The densities of all insect pests were not significantly affected by interaction effect of year × weed management (Table 1).

Population of natural enemies
The effects of year and interaction effect of year × weed management treatment were not significant on densities of all natural enemies (Table 4).The effect of weed management treatment was significant (p ≤ 0.01) on natural enemies densities (Table 4).The highest densities of Coccinella septompunctata, Coccinella undecimpunctata and Orius niger were obtained in CLM treatment (14.4,11.3 and 21.1 no.plant −1 , respectively) and the lowest values (3.3, 2.1 and 5.6 no.plant −1 , respectively) in weed-free treatment (Table 5).The densities of Coccinella    5).

Weed density
The weed species composition in experimental field was presented in Table 6.The grass, broadleaf and total weed densities were not affected significantly by year and interaction effect of year × management treatment (Table 7), but affected significantly (p ≤ 0.01, p ≤ 0.05 and p ≤ 0.05, respectively) by weed management treatment.In weedinfested treatment, the grass weed density was higher than that of broadleaf weed.The grass and broadleaf weed densities decreased significantly at all weed management treatments compared with weed-infested treatment and the highest reductions (76.2 and 89.7%, respectively) were observed in OHW treatment (Table 8).At all living mulch treatments (FLM.VLM and CLM), the grass weed densities were higher than WSM treatment.The broadleaf weed densities in FLM and VLM treatments were higher than CLM and WSM treatments.The total weed density decreased significantly at all weed management treatments compared with weed-infested treatment and the highest reductuion (82.4%) was observed in OHW treatment (Table 8).In CLM, VLM and FLM treatments, the weed densities were not significantly different.The grass, broadleaf and total weed densities in trifluralin treatments were lower than straw mulch and all living mulch treatments (FLM.VLM and CLM).www.nature.com/scientificreports/

Weed biomass
The effects of year and interaction effect of year × management treatment were not significant on grass, broadleaf and total weed biomasses (Table 9).The effect of weed management treatment was significant on grass (p ≤ 0.01), broadleaf (p ≤ 0.01) and total weed (p ≤ 0.01) biomasses.In weed-infested treatment, the broadleaf weed biomass was higher than that of grass weed.At all weed management treatments the grass and broadleaf weed biomasses decreased significantly compared with weed-infested treatment (Table 10).The highest reductions in grass and broadleaf weed biomasses (51.0 and 72.0%, respectively, compared with weed-infested) were observed in OHW treatment.At WSM and all living mulch treatments (FLM, VLM and CLM), the grass weed biomasses were higher than trifluralin treatment.The broadleaf weed biomasses at all living mulch treatments were higher than WSM and trifluralin treatments.The total weed biomass decreased significantly at all weed management treatments compared with weed-infested treatment and the highest reductuion (62.6%) was observed in OHW treatment (Table 10).In CLM, VLM and FLM treatments, the total weed biomasses were not significantly different.The total weed biomasses in trifluralin and WSM treatments were lower than all living mulch treatments (FLM, VLM and CLM).

Sesame seed yield
The effects of year (p ≤ 0.05) and weed management treatment (p ≤ 0.01) were significant on sesame seed yield (Table 1).The sesame seed yield in 2021 (1002 kg ha −1 ) decreased significantly compared with 2020 (1181 kg ha −1 ) (Table 2).The weed-free treatment had the greatest sesame seed yield (1670 kg ha −1 ) and in other treatments, the seed yield decreasd significantly (Table 3).The OHW treatment had the greatest seed yield (1456 kg ha −1 ) among the weed management treatments.There was no significant difference among the seed yields in trifluralin, CLM and WSM treatments.Also, among the living mulch treatments, the seed yields in VLM and FLM treatments, decreased compared with CLM (Table 3).The interaction of year × weed management treatment was not significant on seed yield.

Discussion
The densities of Myzus persicae and Brevicoryne brassicae in 2020 were higher than 2021.This increase could be attributed to higher mean temperatures in 2020 (from April to August; Table 1) compared to 2021.El Fakhouri et al. 32 , also observed that increasing the temperature promoted pea aphid (Acyrthosiphon pisum Harris) population in lentil (Lens culinaris Medikus).Also, the precipitation in 2020 growing season (from May to September) was higher than 2021, that could be another reason for increasing the densities of Myzus persicae and Brevicoryne brassicae.Frank and Liburd 33 also reported that the differences in temperature and rainfall could affect the population densities of aphids.The highest densities of all insect pests were observed in weed-free treatment and reduced significantly in weed-infested and all weed management treatments.These results indicate that the presence of weeds or living mulches in sesame field would reduce the densities of these insect pests.The densities of Myzus persicae and Brevicoryne brassicae in WSM treatment were higher than those in VLM and FLM  34 which reported that the presence of Calendula officinalis enhanced Myzus persicae suppression.Also, in intercropping of cowpea (Vigna unguiculata L. Walp.) and African marigold (Tagetes erecta L.) with Cucurbita pepo, observed that marigold and marigold-cowpea intercropping suppressed the aphids by increase in population of natural enemies 35 .
The lowest densities of Spodoptera exigua and Helicoverpa armigera were also observed in CLM treatment, as Fabrick et al. 36 reported that survival of Lygus hesperus and Bemisia tabaci reduced on French marigold (Tagetes patula L.) plants compared with common bean (Phaseolus vulgaris).However, the biotic mechanisms that confer such repellency are often not well understood 37 and needs more investigation.Moreover, the reductions in insect pests densities in this treatment (CLM) could be related to increase in densities of natural enemies.In WSM treatment, the densities of these insect pests were higher than VLM and FLM living mulch treatments.Frank and Liburd 33 also reported that the densities of whitefly (Bemisia argentifolii Bellows & Perring) in synthetic white and reflective mulches were higher than those of Fagopyrum esculentum and Trifolium repens living mulches.Bruce et al. 30 also observed that in Cucurbita pepo production, the cucumber beetles (Acalymma trivittatum) density in straw mulch treatment, increased compared with those in annual ryegrass (Lolium multiflorum) and white clover (Trifolium repens) living mulches.
The highest densities of all natural enemies (Coccinella septompunctata, Coccinella undecimpunctata and Orius niger) were observed in CLM treatment that caused significant reduction in all insect pest densities.Similar results have been reported in previous studies on pot marigold Calendula officinalis 34 and African marigold Tagetes erecta L. 35,38 .The presence of Calendula officinalis flowers might attract natural enemies into crop fields without any direct effects on natural enemy fitness and improve attack rates on the pest 39 , or it might increase longevity or fecundity of natural enemy 40,41 .Zhao et al. 34 reported that the growth rate and number of Orius sauteri, increased significantly in presence of Calendula officinalis, that consequently increased the Myzus persicae suppression.
All living mulch treatments (FLM, VLM and CLM) had higher densities of natural enemies than WSM treatment.Also, the densities of natural enemies in weed-infested treatment were higher than weed-free treatment.The hypothesis of natural enemy was confirmed by these results, where the increase in densities of natural enemies is the outcome of the increased plant diversity 42 .Frank and Liburd 33 also reproted that the populations of natural enemies in living mulches were higher than control (bare ground) and synthetic mulches.Interplanting of red clover (Trifolium pratense L.) with cucumber (Cucumis sativus L.), increased the densities of natural enemies and decreased the densities of melon aphid (Aphis gossypii) and cucumber beetles (Acalymma trivittatum) 43 .Also, the living mulch of buckwheat (Fagopyrum esculentum) increased the population of natural predator 44 .
At all weed management treatments except FLM, the reductions in grass weed densities compared with weed-infested treatment, were lower than those in broadleaf weed densities.The greatest reductions in grass, broadleaf and total weed densities were observed in OHW treatment that is in agreement with previous studies on Anethum graveolens L. 29 and Glycine max L. 45 .In WSM treatment, the reductions in grass, broadleaf and total weed densities were greater than VLM and FLM treatments that may be attributed to the higher ground cover in this treatment than living mulch treatments (VLM and FLM) 46 .Straw mulches could suppress seedling emergence of weeds through physical barrier created by mulch itself 47 and reduction in light interception by seeds 48,49 .In trifluralin treatment, the reductions in grass, broadleaf and total weed densities were greater than straw (WSM) and living mulches (FLM, VLM and CLM) that could be explained by inhibition effect of trifluralin (as a soil-applied herbicide) on weed seed germination 28 .
At all weed management treatments, the reductions in grass weed biomasses compared with weed-infested treatment were lower than those in broadleaf weed biomasses.The greatest reductions in grass, broadleaf and total weed biomasses were observed in OHW treatment.The Calendula officinalis living mulch (CLM) and wheat straw mulch (WSM) indicated higher weed biomass reduction (weed control efficacy) than Trigonella foenumgraecum (FLM) and Vicia ervilia (VLM) living mulches.In Anethum graveolens, the weed biomass reductions in Trigonella foenum-graecum and Vicia ervilia living mulches were lower than wheat straw mulch 29 .For weed supression, living mulches should have high initial growth rate and competitive ability against weeds 50 .The weed biomass reduction in trifluralin treatment was greater than straw and living much treatments (WSM, FLM, VLM and CLM) that are in agreement with previous findings on Dracocephalum moldavica 28 , Anethum graveolens 29 and Cuminum cyminum 51 .
Among the weed management treatments, the greatest sesame seed yield was obtained in OHW treatment that could be due to the greatest weed biomass reduction (62.6%) in this treatment.Similar results have been reported in previuos studies on Phaseolus vulgaris L. 52 and Anethum graveolens L. 29 , whereas they observed that the grain yields in herbicide application, living mulch and straw mulch treatments, were lower than one-time hand weeding treatment.The highest reductions in seed yield were observed in FLM (44.5% comared to weedfree) and VLM (46.9% comared to weed-free) treatments that may be related to the competition occurred by living mulch treatment as reported in previous studies 29,50 .The sesame seed yield in CLM treatment, enhanced compared with VLM and FLM treatments, while the weed biomasses at all living mulches were not significantly different.These results could be explained by increasing the densities of natural enemies and suppression of insect pests in CLM treatment.The sesame seed yields were not significantly different in WSM and trifluralin application, but the weed biomass reduction in trifluralin treatment (53.1%) was higher than WSM treatment (44.9%).Also the densities of all insect pests in trifluralin treatment were lower than WSM treatment.Therefore, www.nature.com/scientificreports/ the reason for non-significant difference between sesame seed yields in trifluralin and WSM treatments may be related to increase in water content of soil due to straw mulch coverage 29,53 .

Conclusions
Total weed densities and biomasses were not significantly different in CLM, VLM and FLM treatments (living mulches), while the sesame seed yield in CLM treatment was higher than FLM and VLM treatments.The increase in seed yield in CLM treatment compared with other living mulches, could be explained by increase in densities of natural enemies and thereby reduction in densities of all insect pests in this treatment.The densities of all natural enemies in OHW treatment were higher than WSM treatment, which caused reduction in all insect pest densities in this treatment.The highest reduction in total weed biomass and the greatest seed yield were obtained in OHW treatment, threrfore in cropping systems with no labour limitation, the growers could use this method for weed management.The sesame seed yield in trifluralin treatment was not significantly different with CLM and WSM treatments and in sustainable production systems, these non-chemical treatments could be recommended instead of trifluralin for weed management in sesame.One-time hand weeding (OHW) was more efficient for weed suppression than straw and living mulches.This result indicates that developing other weed management methods such as machine powered hoeing and harrowing could improve the weed and pest management, thereby the seed yield in sesame.Also, more investigations are required to investigate the effects of sowing time of living mulch relative to main crop and other types of straw and living mulches in sustainable weed pest management of sesame.

Experimental site
A field study was carried out at East Azarbayjan, Iran (Latitude 38° 53′ N, Longitude 46° 47′ E, Altitude 315 m a.s.l.) during 2020-2021.The experimental area is warm temperate with mean annual precipitation of 375 mm mean annual temperature of 25.8 °C.The 2-year (2020 and 2021) data of monthly total precipitation and mean temperature of the experimental site are presented in Table 11.The properties of soil of experimental field at depth of 0-30 cm are presented in Table 12.

Field practice and experimental procedure
The experiments were arranged as RCBD (randomized complete block design) with eight treatments and four replicates.The weed management treatments were including application of trifluralin (960 g ai ha −1 as preplanting of sesame), using straw mulch of wheat (5.0 t ha −1 ) (WSM), living mulches of fenugreek (Trigonella www.nature.com/scientificreports/foenum-graecum L.) (FLM), bitter vetch (Vicia ervilia L.) (VLM) and calendula (Calendula officinalis L.) (CLM) and one-time hand weeding (OHW) at 50 days after planting (DAP) of sesame.The weed-infested and weedfree treatments during whole growing season also were considered in the experiment (Table 13).In weed-free treatments the weeds were removed manually every day.12 planting rows spaced 50 cm apart with 5 m length (6 m × 5 m) was considered as an experimental plot.
Deep moldboard plowing (25-30 cm) in the spring which was followed by disking before planting was used for seedbed preparation.The soil-applied herbicide, trifluralin (Treflan, EC, 480 g L −1 , Aria Shimi) was applied as pre-planting in to the top soil layer (5-cm); the soil incorporation was done by raking instantly after herbicide use.For trifluralin application a backpack sprayer (Matabi, Goizper Group, Spain) was used equipped with Flood-jet nozzle which calibrated to deliver 265 L ha −1 at 235 kPa.The sowing of sesame seeds was done by hand on 15 May 2020 and 17 May 2021 in soil depth of 2 cm with 18.18 plants m −2 density.The planting row distance was 50 cm and the plant distance over the row was 11 cm.The straw mulch of wheat (5.0 t ha −1 ) was applied between the sesame rows, immediately after planting.The Calendula officinalis, Trigonella foenum-graecum and Vicia ervilia seeds were planted (simultaneously with crop) as living mulches between the sesame rows, at seeding rates of 15, 35 and 40 kg ha −1 , respectively.The growing period of Calendula officinalis, Trigonella foenum-graecum, Vicia ervilia and were 108, 117 and 75 days, respectively.The one-time hand weeding treatment was done 50 DAP of sesame.After planting of sesame and living mulches, the furrow irrigation was done in experimental plots and it was done every 7-day until physiological maturity of sesame.

Ethical approval
This experimental research upon plants complies with relevant institutional, national, and international guidelines and legislation.The seeds of sesame, fenugreek, bitter vetch and calendula were purchased from Tabriz, East Azarbayjan, Iran.

Populations of pests and natural enemies
At 7-day intervals, the densities of pests and natural enemies at different treatments were recorded from their appearance time in flowering stage of sesame (52 DAP; 05 July 2020 and 07 July 2021) until sesame physiological maturity (12 September 2020 and 16 September 2021).To determine the population of aphids in each plot, Myzus persicae and Brevicoryne brassicae were randomly sampled in 20 10-cm twigs in the morning (8:30-11:00 a.m) 32 .For Helicoverpa armigera and Spodoptera exigua at each sampling time, 10 plants were randomly selected in experimental plots and the number of larvae on each plant was counted.
To evaluate the population of natural enemies, Coccinella septompunctata, Coccinella undecimpunctata and Orius niger; the sampling was done every week in 5 × 5 m experimental plots.Coccinella septompunctata and Coccinella undecimpunctata were sampled using a 1 × 1 m quadrat.At each plot, two quadrats were thrown randomly and the ladybird beetles of both species were counted.For data analysis, the mean number of ladybird beetles counted in the two quadrats, were used.A sweep net (38 cm diameter ring and 80 cm light wooden handle) was used for sampling of Orius niger.For a sampling unit, ten 180• sweeps per plot was considered.The insect pests and natural enemies were stored in 70% ethanol in lab condition in vials.

Weed traits
At maturity stage of sesame (121 DAP in 2020 and 123 DAP in 2021) in order to measure grass, broadleaf and total weed density and biomass the weeds were cutted from 1.5 m −2 area (using three 1 m × 0.5 m quadrats).The weed species in the experimental field were included Cyperus esculanthus L., Sorghum halepens L., Setaria viridis L., Portulaca oleracea L., Raphanus raphanistrum L., Convolvulus arvensis L., Silybum marianum L., Cynodon dactylon L., and Solanum nigrum L. The grass, broadleaf and total weed densities and biomasses were measured separately.For weed biomass measurement, the weeds were placed in paper bags; oven-dried at 75 °C for 48 h, and weighed.www.nature.com/scientificreports/

Sesame seed yield
To determine the seed yield, an area of 2 m 2 was harvested from the middle rows of all experimental plots at maturity stage, on 12 September 2020 and 16 September 2021.The mechanical thresher was used to obtain the seed yield (kg ha −1 ) from harvested plants in each plot.

Statistical analysis
A combined analysis of variance (ANOVA) was performed using SPSS software ver.20 based on RCBD with four replications.The data for densities of pests and natural enemies, weed density and biomass and sesame seed yield, met the assumptions of normality and homogeneity of variance and no transformation was needed.For means comparison, the Fisher's Protected Least Significance Difference test (p ≤ 0.05) was used. https://doi.org/10.1038/s41598-023-48978-6

Table 2 .
Densities of insect pests and sesame seed yield in 2020 and 2021.Values given after means are standard errors (± SE).Means within each column with similar letter are not significantly different at p ≤ 0.05 by LSD test.Year Densities

Table 3 .
Densities of insect pests and sesame seed yield at different weed management treatments.The presented data are means for 2 years (2020 and 2021).Values given after means are standard errors (± SE).Means within each column with similar letters are not significantly different at p ≤ 0.05 by LSD test.aTheabbreviations of the weed management treatments were presented in Table13.

Table 4 .
Analysis of variance for effect of weed management treatments on densities of natural enemies.ns, * and **: non-significant and significant at p ≤ 0.05 and p ≤ 0.01, respectively.

Table 5 .
Densities of natural enemies in different weed management treatments.The presented data are means for 2 years (2020 and 2021).Values given after means are standard errors (± SE).Means within each column with similar letters are not significantly different at p ≤ 0.05 by LSD test.aTheabbreviations of the weed management treatments were presented in Table13.

Table 6 .
Common name, scientific name, family name and morphology of identified weed species in experimental field.

Table 7 .
Analysis of variance for effect of weed management treatments on grass, broadleaf and total weed density.ns, * and **: non-significant and significant at p ≤ 0.05 and p ≤ 0.01, respectively.

Table 8 .
Effect of different weed management treatments on grass, broadleaf and total weed density.The presented data are means for 2 years (2020 and 2021).Values given after means are standard errors (± SE).Means within each column with similar letters are not significantly different at p ≤ 0.05 by LSD test.aTheabbreviations of the weed management treatments were presented in Table13.Vol.:(0123456789) Scientific Reports | (2023) 13:21586 | https://doi.org/10.1038/s41598-023-48978-6

Table 9 .
Analysis of variance for effect of weed management treatments on grass, broadleaf and total weed biomass.ns, * and **: non-significant and significant at p ≤ 0.05 and p ≤ 0.01, respectively.

Table 10 .
11fect , that are in consistent with findings of Frank and Liburd33.They found that in Zucchini, (Cucurbita pepo L.), the aphid density in synthetic white mulch treatment increased compared with living mulches of white clover (Trifolium repens L.) and buckwheat (Fagopyrum esculentum Moench).Conversely, in bush bean (Phaseolus vulgaris), Gill et al.11observed that aphid densities in living and straw mulches and unmulched control, were not significantly different.The lowest densities of Myzus persicae and Brevicoryne brassicae were observed in CLM treatment that is in agreement with Zhao et al.
of different weed management treatments on grass, broadleaf and total weed biomass.The presented data are means for 2 years (2020 and 2021).Values given after means are standard errors (± SE).Means within each column with similar letters are not significantly different at p ≤ 0.05 by LSD test.aTheabbreviations of the weed management treatments were presented in Table13.treatments

Table 12 .
Physico-chemical properties of the soil of experimental area in depth of 0-30 cm.